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Abstract 

We study the spin-dependent transport through a quantum dot coupled to two fer- 
romagnetic electrodes using the equation of motion method for the nonequilibrium 
Green's functions. Our results show that the conductance and the density of states 
(DOS) are strongly dependent on the configurations of the magnetic electrodes. 
In parallel configuration of magnetic electrodes the conductance is affected by the 
spin-flip process and the Coulomb repulsion on the dot. The Kondo peak for spin- 
dependent transport is splitted into two peaks by the spin-flip process. The locations 
of the two peaks are symmetric about no spin-flip peak and the separation of the 
splitting is dependent on the strength of the spin-flip parameter R. This effect may 
be useful to realize the spin-filter device. 
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1 Introduction 



In metallic bulk the dilute magnetic impurities are screened by the itiner- 
ant electron spins and cause anomalous resonant scattering of conducting 
electronsfl]. It shows up a minimum in electric resistivity at the low tem- 
perature. In the past decade Kondo phenomena in quantum dots coupled to 
metallic leads has attracted much attention[2,3,4,5,6,7,8] and has been studied 
intensively in both theories and experiments. There is a fundamental difference 
between the Kondo effect in quantum dots and the Kondo effect in alloys with 
magnetic impurities. In the latter case the dressed localized magnetic moment 
increases the scattering cross section and as a result there is a minimum in the 
resistivity, whereas for the former case the resistivity depends on the parity 
of the level occupied by electrons. However, the physics responsible for these 
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resistances is similar. In both cases the effect is due to the spin correlated 
transport at the Fermi level. 

The Kondo effect has been understood in the structure that a quantum dot 
coupled to the normal electrodes in equilibrium and out of equilibrium[9,10,ll,12,13,14,15,16,17,18]. 
It is a consequence of high order perturbation and co-tunnelling process caused 
by the electron-spin interaction between electrons in the quantum dot and in 
lead. The Kondo effect through quantum dots coupled to external leads can 
even induce current in the Coulomb blockade regime by making spin cor- 
related singlet and triplet with electron spins outside the quantum dot [3,4]. 
The conductance is related to the number of electrons in the dot. Usually the 
Kondo effect appears when the dot has odd number of electrons (total spin 
1/2), but is absent when the dot has even number of electrons (total spin 0). 
The conductance or the density of states (DOS) of the dot is enhanced by the 
Kondo effect at low temperature. However, for the ferromagnetic-quantum 
dot-ferromagnetic system, the question is open. We want to know what is 
the manifestations of the Kondo effect when the magnetic moments of the 
electrodes are taken into account in finite and infinite Coulomb repulsion. 

In this work, we investigate the spin-dependent transport in an interacting 
quantum dot coupled to two ferromagnetic electrodes. Different from the con- 
ventional Kondo problem in normal electrodes, i.e., quantum dot-normal elec- 
trodes system, here we deal with the strong electronic correlation, which is 
sensitive to the relative orientations of magnetization between the two elec- 
trodes. The paper is organized as follows. In Sec. 2 we describe the model of 
the system and then calculate the nonequilibrium Green's functions using the 
equation of motion method. The general expression for the electric current 
flowing through the dot in terms of nonequilibrium Green's functions is given. 
In Sec. 3 we present numerical results with and without the spin-flip process 
in the dot for the arbitrary Coulomb repulsion U. We have some conclusions 
in Sec. 4. 



2 Model 

The Hamiltonian for the quantum dot coupled to two ferromagnetic electrodes 
reads 

H= £kaaat acT a ka(T + £ddtd a + U df t^rffc^ + R{d~\d l + dj"d T ) 

k;a£L,R;a ° 

+ E (V k «at aa d a + h.c). (1) 

k;a£L,R;cr 
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Here the single-particle energy Ed in quantum dot is double degenerate in the 
spin index a (<7=f|, ±), and the Coulomb repulsion interaction in the dot 
is U . The first term in Hamiltonian describes the free electron in magnetic 
electrodes, the second and third terms represent the correlated level of the 
quantum dot, the fourth term is to describe spin-orbit coupling which may 
cause the spin rotation of an electron in the quantum dot, and the last term 
is the spin-dependent hybridization of the quantum dot with the magnetic 
electrodes. 



Since the spin quantization axes in the electrodes are fixed by the internal 
magnetization of the magnets, an electron is in a superposition of spin-up 
and spin-down states as it tunnels into and out the dot. In calculations we 
have to taken into account the coherent processes. Technically, we take a 
tr ansf ormat ion [1 8 , 1 9] , 

d a = -^=(c a - crc„). (2) 



In terms of Eq. (2) the Hamiltonian becomes 



k;a£L,R;cr & 

+ 2 (~lKV ka at aa c a - -^=Vk aa 4 aa c<r + h.c), (3) 

where e ca =Ed + <jR for spin-cr. The current through the left electrode can be 
calculated in terms of nonequilibrium Green's functions with the tunnelling 
matrix 

V ka = ^=\ ^ Vkai I • (4) 



V2 



The current is 
2e 



J = i^Y, j ' Tr[Vk L G< kL {u)], (5) 



where [G^ L (t)} a(T '=i(a kLa .(0), c a >(t)) is the less Green's function of the left elec- 
trode. 

Now we apply Dyson's equation to express the current J by the Green's func- 
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tions G of the dot, 

J = \ I ^^{ ri H G< H + Me)[G R (u>) - G A (u)]}, (6) 

where the less Green's function of the dot [G < (t)] aa '=i(c^,(t), c CT (0)). The re- 
tarded (advanced) Green's function of the dot [G R ^ A \t)] aa i=^f i({c CT (£), cj(0)}), 
respectively. The upper sign is for retarded Green's function and the lower sign 
is for advanced Green's function. / Q (e) is the Fermi distribution function in 
the a electrode, and 



r>) = i ( r?H + rf M ry( w ) - r?M 
2 ^ r^(cj) - t«(lu) t^(u) + r^) 



(7) 



is the linewidth matrix with T"(u)=2nJ2k I Vkaa- | 2 ^(^ — £kaa)- The spin 
dependence of T"(e) originates from the bulk magnetization of the electrodes. 

In order to determine the current J, we should first determine the Green's 
function of the dot. As we know the retarded Green's function G R is a con- 
jugation of the advanced Green's function G A , we need only to calculate G R . 
With the notation as usual G R a , =< c CT ,c+, >= -i J °° e iwt ({c a (t), c+,(0)})<Zi, 
we apply the equation of motion method again to calculate the Green's func- 
tion G R a ,. The first equation of motion for the retarded Green's function G R a , 
is 



< c CT , c+ >= ({c CT , c+ })+ < [c CT , if], c+ > . (8) 



Using the Eq. (8), G R a , can be written as 

(w - e«,)G&, = <W + tfc© + £ -^=n; CT r^' - £ ^n;,r-', (9) 



where the second order function G^, , T^' and T^' are G^j, c a c£ca, c£, ^> 
> r fca =< c+, >, and Tf^' =« a fcQ(f , c+, >, respectively. 

There are three new Green's functions on the right hand side of Eq.(9). The 
equation of motion for T^' and are close since the only new Green's 
functions are G R a , and G R a ,, 

Tfca = -^ V kaag r la G R a , - -j=V kaa g r la G R a , , (10) 
Tfca = ^^toa&^ - -j^VkaaQlaG^, , (11) 
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where g{ ff =(u - e kaa + io + ) 1 and g r 2 ^{u - e koi& + io + ) 1 . Note that g r la =g r 2a . 

(2) 

However, the equation of motion for the new Green's function G^,, which 
is generated by the Coulomb interaction and involves a two-particle Green's 
function, is not close. We have to truncate the equations. An approximate solu- 
tion, valid for temperatures higher than the Kondo temperature [20], ksT^ — 
(f/r)( 1//2 Wp(— 7i\fi — e a |/r), is obtained by neglecting terms in the equation 
of motion for G^, which involves correlations in the leads[10]. We can write 
out G^l, as follows, 



k;a V 1 k;a V 1 k;a V 1 

-£4^uf, (12) 

k;a V2 

with the new notations: rf=<^ a ka actca, c+, >, 1^=^; a ka5 c^Ca,c^, >, 

Tg 2 ^^ Co-c^afcao-, c+, ^>. n„=(c^Cv) must be calculated self-consistently[ll]. 

As for r| 2 \ we adopt the truncation as Y. Meir et al[9,10] and write out vf^ 
explicitly as, 

= —^V kacT gl a G^, — -j=V kaa g[ a f(e kaa )G^ a ,, (13) 
r 2 2) = --^j^kaaglaG^l, + -j= V ka ^g r 2 ^f (e kad )G^ a , , (14) 



1 0~ 1 

^3 — — -^VkaadloG aa ' ~ ~^J^k(xa9w^ W + ^J^kaadla) '( £ kaa)G ' aa , 

+^V k * a ,gZJ(e ka ,)Gl„ (15) 



p(2) _ °" T/ » r ^-,(2) [ Tr* r .^,(2) ' T t^* r r/ 

-^V k * aa g r 4 J(e kaa )Gl„ (16) 

— ^^fca^a^W ~~ -^Vka&glaJ " ( £ kaa) G^ a , , (17) 
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F (2) _ 
1 6 — 



(18) 



with the notations as g^=(u — e ca — e cd - + e ka ^ — U + io + ) 1 , g\ a ={oj — e ca — 
tea + £kaa -U + io+y 1 , g^=(uj - e ca + s c5 - e ka5 + io + y l , and g r %a ={uj - 
£ C a + £ca ~ £kaa + io + )~ l . Notice that g\ a =gl a . 

Defining E i(7 =Efc ; a \ I V kaa \ 2 g r ia and E{ T =Efc ;a \ I ^ P f(^kaa)g r ia , we have 
G^ 2 j, in the form of 



(u - e ca - U)G®, = n„5 aa > + (S 1(J + E 2s + E 35 + S 4ct + S 5S 
+S 6CT )GS + (<rE 3s + - (E^ + Ei + 



(19) 



Making the truncation as above, the equation for G a ^, is close. It now only 
involves the Green's functions G*®,, G^,, G R a , and G§ a ,. One can solve them 



self-consistently. With the notations of e C(T = 



(e 

C C(T 


0^ 






') 


I 


£ ca j 









Ei| + Yi2i Eij_ — E21 
Sij — E21 Eif + E2f 



E^ + Eg| s{| E^ ^ 



Ei-f + E24 + S51 + 

Eif + + Esf + E6j, 



, and 



<f). 



(y{f) , 

Zj 3 | -|- Zj 4 j_ Zjgj Zj 4 | 



V (/) V (/) V (/) , y(f) 



, defining E 4 =wl — e ccr — C/l — E 2 — E 3 , one 

can write the Green's function G R in a compact form as below 

G R = [ul - e ca - E + Ut^iti + Ex)]-^! + Urvt A v ). (20) 



In the infinite U limit, we recover the result in Refs. [9,18], 
G R = (col - e cc - E - E^-^i - n). 



(21) 



To determine the lesser Green's function of the dot G < , we use the Ng's 
ansatz[21] in our case. The interaction lesser and greater self-energies are as- 
sumed to be of the form E < ' > =Eq A, where A is a matrix to be determined 
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by the condition of S < -S > =S i? -S A . This ansatz is exact in the noninteract- 
ing limit (£7=0) and guarantees automatically the current conservation law. 
As a result one obtains S < =Sq [S^ — Sq] -1 ^^ — S A ]. Using this results, G < 
can be obtained by Keldysh equation G < =G R T I < G A . Substituting the ex- 
pressions of the Green's functions of quantum dot into Eq. (6) and defining 
T I < =T R (T L + r^) -1 ^^ — £ A ), one obtains an expression of the tunnelling 
current 

J = I / ^ T ^ L G R Z<G A ][f L (e) - f R (e)\. (22) 

This expression generalizes the well know current formula to the spin-dependent 
Anderson model with additional spin-flip relaxation and allows one to describe 
the coherent spin transport through an interacting quantum dot coupled to 
magnetic electrodes. 



3 Numerical results 

In this section, we will present the numerical results of spin-dependent DOS 
and conductance. For simplicity, we neglect the energy dependence of the tun- 
nelling matrix elements and consider parallel and antiparallel configurations. 
If the magnetic moments in the electrodes are parallel, the spin-majority elec- 
trons are assumed to be spin-up and the spin-minority electrons are spin-down. 
In the antiparallel configuration, the magnetization of the left electrode is up 
and the magnetization of the right electrode is down. With these assump- 
tions the coherent spin transport parameters can be conveniently expressed 
by introducing magnetic polarization factors pi and p^ for the left and right 
barriers, respectively. Therefore, r^=r (l ±p L ) and T R ^=aT (l ±Pr) are 
for the parallel configuration, and rwjN=r (l ±Pl), r^^oT^l TPr) for the 
antiparallel configuration. To describes the coupling between the quantum dot 
and the left electrode without internal magnetization. We express all parame- 
ters in the unit of r in the following calculation, a denotes tunnel asymmetric 
factor of the left and right barriers. In this work, we assume the symmetric 
barriers, i.e., a—1, Pl=Pr=P- 

The spin-resolved DOS in the dot are calculated via the relation p-m) = 
— l/nIm{G R ^fG^^fG R +Gf^}. The Kondo resonance for each spin is clearly 
manifested by the peak at oo =0 in the DOS for both parallel and antiparallel 
configurations. However, the shape of the peak is sensitive to the magnetic 
configurations of the electrodes. As observed from Fig. 1(a), in the parallel 
configuration the spin-up DOS (solid line in Fig. 1(a)) is dramatically differ- 
ent from that of spin-down (dotted line in Fig. 1(a)) in the infinite U limit, 
while in the antiparallel configuration the DOS for spin-up and spin-down are 
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CO CO 

Fig. 1. Density of states in the parallel (a) and antiparallel (b) configurations in the 
infinite U limit, the solid line is for spin-up and the dotted line is for spin-down. 
Here e d =-2, X S T=0.001, p=0.5, R=0. 

identical, see Fig. 1(b). In Fig. 1(b) the Kondo peak for spin-up and spin-down 
is not affected by the magnitude of magnetic polarizations p. As we know that 
the spin current of electrons through the antiparallel configuration cannot 
be polarized, whereas can be polarized through the parallel configuration [22]. 
Therefore, the DOS for spin- up is different from that of spin-down in the par- 
allel configuration and is the same as that of spin-down in the antiparallel 
configuration. 

1.0 

0.8 

CO 0.6 
O 
Q 

0.4 
0.2 
0.0 

-10 -8 -6 -4 -2 2 4 6 8 10 

CO 

Fig. 2. Spin- up density of states in the antiparallel configuration for [7=5 (solid line), 
6(dotted line), and 7(dashed line). Here R=0, the other parameters are the same as 
in Fig.l. 

In the limited Coulomb repulsion U, the spin-up DOS in the antiparallel con- 
figuration are plotted in Fig. 2. It can be seen from Fig. 2, there are two 
broad maxima in DOS, which are centered at the positions of the dot level Ed 
and Ed + U, respectively. However, for infinite U limit shown in Fig. 

1, there is only one broad maximum. It is consistent with physical intuition 
since for the infinite U there is only one level in the dot but are two levels for 
the finite U. This result agrees with other authors' work[18,12]. The Kondo 
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Fig. 3. Total density of states in the antiparallel configuration with v=0 (solid line), 
0.5 (dotted line), and f (dashed line), respectively. Here C/=6, the other parameters 
are the same as in Fig.f. 

peak in DOS is exactly located at the Fermi level, see the sharp peak in Fig. 
2. This is the equilibrium situation (the applied bias is set to be zero). When 
the voltage is applied (the left chemical potential is different from that of the 
right), the Kondo peak splits into two peaks with lower intensities at positions 
of the left and the right chemical potentials. This behavior is shown in Fig. 
3 for three different values of the bias voltages. As the bias voltage increases 
the splitting of the Kondo peak becomes wider and is equal to the difference 
of the bias. The Fig. 3 is plotted in a narrow energy range around the Kondo 
peak, the inset of the figure gives a wider view in energy. It is evident that 
the peaks occur at the Fermi levels of the two electrodes. The solid line is for 
zero bias, there is only one peak at uj = 0. The dotted line is for bias v=0.5 
and the second peak appears at the position uj = 0.5. The dashed line is for 
v=l, the splitting peak is at uj — 1. This behavior is in accordance with other 
results[23,24,25]. 

The interesting feature in this work is that there exists a splitting due to the 
spin-flip transition processes. In the preceding discussions there is always a 
Kondo peak at the energy uj = 0. However, if we set spin-flip parameter R 
not to be zero, as shown in Fig. 4 for the DOS of spin-up in the antiparallel 
configuration, the original Kondo peak splits into two peaks. The two peaks 
appear at the positions of uj = ±R, respectively, and the original peak at uj = 
disappears. In Fig. 4 the solid line is plotted without spin-flip transition, the 
dotted line is for spin-flip transition parameter R = 0.2 and the dashed line 
is R = 0.5. The inset in Fig. 4 is also a global view of the Kondo effect. 
The splitting is due to the coupling of spin-flip transition processes between 
the dot and the magnetic electrodes and can be expressed by the high order 
self-energy. 
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Fig. 4. Spin-up density of states in the antiparallel configuration with R=0 (solid 
line), 0.2 (dotted line), and 0.5 (dashed line), respectively. Here U=6, the other 
parameters are the same as in Fig.l. 
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Fig. 5. The conductance in the antiparallel and antiparallel configurations, (a) in the 
antiparallel configuration for the infinite U limit and different temperatures. The 
spin-flip parameter R=0. (b) in the antiparallel configuration for different R, here 
U=6 and KbT=0.001. (c) in the parallel configuration for different R , the other 
parameters are the same as in (b). (d) the conductance in the parallel configuration 
for the spin-up and spin-down, respectively, KbT=0.001, the other parameters are 
the same as in (a). 

In Fig. 5, we show the linear response conductance as a function of energy 
of the dot in the parallel and antiparallel configurations. In the antiparallel 
configuration the conductance is similar to the normal case, i.e., the Kondo 
effect broadens the conductance peak, see Fig. 5(a). The peaks are similar 
in shape for different temperatures. They only shift a little to each other. 
As temperature increases the peak shifts to the higher energy. This is caused 
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by the shift of the Fermi energy for different temperatures. Due to the DOS 
for spin-up and spin-down are the same in the antiparallel configuration, as 
seen in Fig. 1, the conductances for spin- up and spin-down are overlapped 
in this configuration. However, when U is finite, for instance £7=6, the con- 
ductance of spin-up differs to that of spin-down, see solid line in Fig. 5(b). 
If the spin-flip process on the quantum dot is involved, this shift becomes 
more dramatically, see the dotted line and dashed line in Fig. 5(b). For a large 
value of R the amplitude of conductance is small. In the parallel configuration 
the conductances for spin- up and spin-down are different [22] , see Fig. 5(d), 
the resonance of spin-up conductance (solid line) has an apparent shift to the 
spin-down (dotted line) due to different densities of states for them [see Fig. 
1(a)]. The existence of the spin-flip transition R reduces the spin-down con- 
ductance under the finite U condition, see Fig. 5(c). Increasing the value of 
R the spin-up conductance becomes higher and the spin-down conductance 
becomes lower. This effect is caused by the correlation between the dot and 
the ferromagnetic electrodes. For a large enough R the peak of conductance 
of spin-down is suppressed and the spin-up peak is enhanced. We expect this 
result may be important in the exploiting the spin-filter devices. 



4 Conclusion 



Using the Anderson model, we study the Kondo effect in the ferromagnetic- 
quantum dot-ferromagnetic system. We find that the spin-dependent transport 
in such a system is sensitive to the alignment of the magnetic moments in 
electrodes. Due to the external applied bias or the spin-flip process on the 
dot, the original Kondo peak splits into two peaks. The conductances in the 
parallel and antiparallel magnetic configurations are affected by the spin-flip 
transition and the Coulomb repulsion on the dot. We expect that the effects 
can be useful to realize the high spin polarization devices. 
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